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Abstract 

We consider a dynamical brane world in a six-dimensional space- 
time containing a singularity. Using the Israel conditions we study the 
motion of a 4-brane embedded in this setup. We analyse the brane 
behaviour when its position is perturbed about a fixed point and solve 
the full nonlinear dynamics in the several possible scenarios. We also 
investigate the possible gravitational shortcuts and calculate the delay 
between graviton and photon signals and the ratio of the correspond- 
ing subtended horizons. 
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1 Introduction 



It has recently been argued that there might exist some extra spatial di- 
mensions, not in the traditional Kaluza-Klein (KK) scheme [1] where extra 
dimensions are compactified on a radius of the order of Planck scale, but in 
a scenario where they could be large. Certainly, the extra dimensions idea is 
not completely new but it is also related to string theories. Among the exist- 
ing models the ten-dimensional E 8 x E 8 heterotic string theory seems to be 
the most acceptable candidate to describe our world. This theory is related 
to an 11-dimensional spacetime, where the 11th dimension is compactified 
via a Z 2 orbifold symmetry. In this setup the standard model particles are 
confined to a four-dimensional spacetime, while the gravitons propagate in 
full spacetime. 

Large extra dimensions were also recalled to solve the hierarchy problem 
[2]. However, it was the work of Randall and Sundrum [3, 4] which suscitated 
a renewed interest on these grounds. In this model our world is identified 
with a domain wall in five-dimensional anti de Sitter (AdS) spacetime. In 
their first paper, Randall and Sundrum proposed a mechanism to solve the 
hierarchy problem by a small extra dimension, while in their second paper, 
the braneworld with a positive tension was investigated. 

In this paper we consider a six-dimensional model where the bulk contains 
a singularity and is bounded by a four-dimensional brane containing matter. 
This brane is a thick brane with the usual three infinite spatial directions 
and an additional coordinate compactified in a small radius that can be of 
Planck size like in the KK models [1] or even of submilimetric size as in 
the Arkani-Hamed, Dimopoulos and Dvali (ADD) models [2], which makes 
it unavailable to the brane observers. In this way, this "effective" picture 
corresponds to our universe. The static case has been already studied in 
[6] where AdS-Schwarzschild and AdS-Reissner- Nordstrom embeddings were 
considered. The case of domain walls moving in bulks containing non-charged 
singularities has also been studied in [7, 8]. Brane models in AdS space 
with Schwarzschild singularities have been used to understand the AdS/CFT 
correspondence and look like promising theoretical models [9]. 

The motivation for a six-dimensional model comes from the fact that 
spacetimes with more than one extra dimension can allow for solutions with 
more appealing features, particularly in spacetimes where the curvature of the 
internal space is non-zero. More extra dimensions also relax the fine-tunnings 
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of the fundamental parameters. On the other hand, the ADD model [2], 
which relates the Planck mass to the fundamental mass in (4 + ra) dimensions 
through the size of the transverse dimensions, foresees the existence of n > 2 
extra dimensions. The possibility n — 1 remains excluded from experimental 
bounds since there are no observed violations of the Newton's law at distances 
of the order of the solar system. 

According to the model, photons are confined to the brane and gravitons 
can propagate in all five spatial directions, then there is a possibility that 
gravitational fields while propagating out of the brane speed up due to the 
warped bulk geometry reaching farther distances in smaller time than light 
propagating inside the brane, a scenario that for a resident of the brane (as 
ourselves) implies shortcuts [5, 6, 10, 11, 12]. 

The existence of shortcuts implies that causality is violated from the 
point of view of a four- dimensional observer; however, it is perfectly defined 
through iV-dimensional geodesies for an observer in the bulk. Thus, two 
points apparently causally disconnected in four dimensions could be causally 
connected by gravitational shortcuts. This fact accounts for a possible solu- 
tion of the horizon problem. As it was pointed out in [13], if high redshifts 
were available, shortcuts appearing before nucleosynthesis are serious medi- 
ators of homogenization of the matter on the brane and then they could pro- 
vide an explanation alternative to inflation to this problem. Other alternative 
approaches have also considered a varying speed of light [14]. Moreover, if 
inflation took part on the brane, the causal structure is definitely changed by 
shortcuts possibly leading to a non-usual period of causal evolution of scales, 
what could be responsible for distinct predictions in the cosmic microwave 
background structure for inflationary models. 

The paper is organized as follows. In section 2 we describe the six- 
dimensional model considered along this work, derive the equation system 
governing the brane motion from the point of view of a bulk observer and 
write the geodesic equation corresponding to the shortest graviton path in 
the bulk as well as the expression for the delay between graviton and pho- 
ton flight times. Section 3 is devoted to studying the brane behaviour when 
its position is perturbed about a fixed point. In section 4 we provide the 
solutions of geodesic equation and brane equation of motion in the different 
scenarios resulting from the combination of the parameters appearing in the 
model and show the existence of several shortcuts. We also find the time de- 
lay and the ratio between the horizons subtended by gravitons and photons. 



3 



Finally, we discuss our results in section 5. 



2 The Brane Cosmological Model 

We consider a six-dimensional model described by the following metric 

ds 2 = —n 2 (t, y, z)dt 2 + a 2 (t, y, z)dY? k + b 2 (t, y, z)dy 2 + c 2 (t, y, z)dz 2 , (1) 

where dT? k represents the metric of the three-dimensional spatial sections 
with k — — 1, 0, 1 corresponding to a hyperbolic, a flat and an elliptic space, 
respectively 

The matter content on the brane is directly related to the jump of the 
extrinsic curvature tensor across the brane [8, 15]. This relation has been 
derived in the case of a static brane in a previous work [6] . Here we generalize 
our result for the Israel conditions to include the case of a brane moving 
with respect to the coordinate system. Its position at any bulk time t will 
be denoted by 

z = K(t) . (2) 
The extrinsic curvature tensor on the brane is given by 

K M n = Vm Vl n N , (3) 
where h A is a unit vector field normal to the brane worldsheet and 

Vmn = 9mn - n M n N (4) 

is the induced metric on the brane. 

In order to compute the components of n A , we use the relations 

g MN h M h N = l , g MN n M u N = 0, (5) 

where we have introduced the unit velocity vector corresponding to the brane, 
which reads 

The relation between dt (the bulk time) and dr (the brane time) can be 
found from the induced metric, 



d? 2 

u,J induced 



- n 2 (t, K{t)) - c 2 (t, n{t))z 2 \ dt 2 + a 2 (t, K(t))dE 2 w 
-dr 2 + a 2 (r)cffij 4) , (7) 
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where a dot means derivative with respect to the bulk time t. We obtain 
dr = n(t,H(t)) 



Thus, (6) can be written as 

u A = ^{l, 0,0, 0,0, ft} , (9) 
and from (5) we can easily obtain 

n A =\ CU 2 , ,0,0,0,0, 1 = , } ^ (10) 

Now we can calculate the components of the extrinsic curvature tensor 
by substituting (10) into (3). The non-zero components are 

" n 2 \ L n*(t,K(t))) 

x\n + — - ft - - 2- - ft 2 2 ft 3 — , 11 

[ c z \n cj \ n c J n A J 

K« = KK°, K° = ~Kk», (12) 

c V n 2 (t,ft(t)) / \ a an 2 J 1 ' v ; 

= l^- ^^ V^ + ^ft (14) 
5 n 2 (t,ft(t))J \6 + 6n 2 J ' 1 j 

*S = (15) 
where all the coefficients take values on the brane. 



2.1 The Israel Conditions 

The energy-momentum tensor on the brane located at z can be written as 
t mn = — — — {(p + p)umu n + PVmn} ■ (16) 
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We also define a tensor T AB as 



f AB = T AB - l -T VAB . (17) 



The Israel junction conditions [16] are given by 



[K^\ = -4)7^ , (18) 



where the brackets stand for the jump across the brane and = e^e„K AB , 
where e„ form a basis of the vector space tangent to the brane worldvolume. 



The non-zero components of (17) are given by 



rpO 
1 


= _± 

C 


[3p + 4p] 


(19) 


rp6 


= 7^f ° 




(20) 


n 


= 

4c 


> J 5 — 1 i > 


(21) 


rh<6 

1 6 


c 2 
n 2 




(22) 



The left-hand side of (18) can be calculated taking into account the mirror 
symmetry across the brane 

[K, u ] = K^HifY) - K^Tlif)-) = -ZK^Kit)) . (23) 

At this point it is convenient to choose a specific bulk metric of the form 
(1) satisfying six-dimensional Einstein equations. This is given by 



where 



ds 2 = -h(z)dt 2 + a 2 (z)dT, 2 k + h~ 1 (z)dz 2 , (24) 



a(z) = Z - , (25) 



,2 dr ' 2 



and 



d ^ = T^yT2 + r d %) + (! - kr W , (26) 

z 2 M O 2 
M*) = * + |f-f + (27) 
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with l~ 2 oc —A (A being the cosmological constant, which can be positive or 
negative) and M and Q 2 are constants. 

We should stress that dy can be written as dy = R c dip, with if an angle 
with the usual periodicity 2n. Thus, the coordinate y is compactified under 
some mechanism such that its radius of compactification R c is small enough 
to evade experimental detection [1, 2]. This makes the local observers have 
the picture of living on an effective three-dimensional brane. 

Metric (24) contains a singularity located at z — 0. It is valid on the 
z < TZ(t) parts of surfaces of constant t and its reflection, by the Z 2 orbifold 
symmetry, is valid on the z > 1Z{t) parts. If M = and Q 2 = 0, then (24) 
is simply the metric of de Sitter or anti de Sitter spacetime according to the 
sign of I 2 . 

With this Ansatz the Israel conditions (18) reduce to only two equations, 
which read as 



where again all the metric coefficients must be evaluated on the brane. Sys- 
tem (28) describes the full nonlinear dynamics of the brane embedded in the 
static bulk (24). 

2.2 The Geodesic Equation and the Time Delay 

We consider two points on the brane ta and tb- In general there are more 
than one null geodesic connecting these points in the 1 + 5 spacetime. The 
trajectories of photons must be on the brane and those of gravitons may be 
outside. The graviton path is defined equating (24) to zero. Since we are 
looking for a path that minimizes t when the final point r# is on the brane, 
the problem reduces to an Euler-Lagrange problem [5]. Then as in [17] the 
shortest graviton path is given by 




(28) 



a' TZ a 
a h 2 a 





. 



(29) 
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We can also calculate the time delay between the photon travelling on 
the brane and the gravitons travelling in the bulk. Since both signals cover 
the same distance [13], 



Jo 



t f dt n 



ft(r 7 ) Jo n g {t g )\ 



h{n a 



K g {t) 2 

h{n g ) ' 



(30) 



the difference between photon and graviton flight times measured by an ob- 
server on the brane can approximately be written in terms of the bulk time 
t as [17] 



At 



~ ll(t f ) f 
Jo 



dt 



TZ g (t)\ 



{ g} h{K g ) K(t) dt 



(31) 



It is also interesting to look at the ratio between the horizons subtended 
by the photons travelling on the brane and the gravitons travelling in the 
bulk. This ratio uses the same quantities previously quoted for the time 
delay, 
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Ml! 

h(K g ) 



lo 



*/ dt dr 
Tl(t) dt 



(32) 



3 Brane Fluctuations 



In this section we study perturbatively the brane behaviour at a fixed point. 
We define the background solution as the case where the brane position is 
frozen, i.e. 

z = Tl{t) = TZ = const. (33) 
Then the system (28) reduces to 



Yh 

a 



K 



(6) 



(3p + 4p) 



"(e) 



p. 



(34) 
(35) 
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The equations of motion for the brane fluctuations can be obtained lin- 
earizing the exact equations of motion (28) about the "equilibrium" config- 
uration. They read as 



57Z 1 „ (h'\ k 2 



I ■ i u j> r ! uji) 

(36) 



a 



From the explicit expressions of h (with Q — 0) and a we can find 

S (7) = -^«K = myw (37) 

5* (9 = {-lk (zp+ip)2 + k l y % } ro 5 mlsn ■ (38) 

As we can see from (36), the brane fluctuations are bound to the matter fluc- 
tuations. However, when we have the case of adiabatic matter fluctuations 

5p = v 2 5p, (39) 

we can derive an equation for the brane fluctuations with an appropriate 
linear combination of (36) 

+ (mi + (3 + 4 v!)mi)5n = 0, (40) 
while the matter fluctuations will be related to the brane ones by 

^5p = ml5n. (41) 
Assuming an equation of state p = up, (40) finally reads 

SU + ! T^A - (7 + M(3 + ) m = . ( 12) 



h 2 U + f -# V " ' 64 
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The constant M can be expressed in terms of 1Z and uj from the Israel con- 
ditions for a static brane [6] as 

M = irk{^ + Z)k+ ^ + l) ^}- (43) 

In the case of a flat domain wall (k = 0, uj = —1), the term inside braces 
in (42) can be interpreted as an effective cosmological constant on the brane, 



4 

V = -A e// , (44) 



1 K (6) _ 2 _ 



n 2 64 

which vanishes when we take into account (35). This corresponds to the 
Randall-Sundrum condition [3, 4] in six dimensions. Then, the effective cos- 
mological constant on the brane is zero. 

Therefore, we are left with the equation of motion for a free scalar field as 
in five dimensions [15, 18]. In this case the brane is at rest and the equilibrium 
position 71 can be chosen arbitrarily. However, a small departure from this 
position results in a runaway behaviour. 

This result is compatible with the numerical solution of the full nonlinear 
system (28) as we will see in the next section. 



4 Nonlinear Brane Dynamics 

The system (28) describing the brane dynamics can be numerically solved 
for several combinations of M, Q 2 , k and I 2 . When M, Q 2 and k vanish, 
the solution for 1Z(t) is a constant or a linear function in t depending on 
the given initial condition for TZ, and we verify the behaviour found through 
brane fluctuations. 

For M and Q 2 non-zero we have solved (28) in the typical cases of a 
domain wall (uj = —1), matter (uj = 0) and radiation (uj = 1/3) dominated 
branes. In figure 1 we show all the possibilities for the bulk metric according 
to the form that h(z) takes due to M, Q 2 , k and I 2 combinations, which 
we study in this section. Some of our results are also illustrated in what 
follows together with the solution of the geodesic equation (29) in order to 
verify the possibility of having shortcuts. We have also calculated the time 
delays and the ratio between graviton and photon horizons for the examples 



10 




(g) (h) (i) 



Figure 1: h(z) for all the possible combinations of M, Q 2 and I 2 , including 
k = 1 (dashed lines), k = (solid lines) and k = —1 (dotted lines) cases. 



of shortcuts appearing in this paper, these are shown in Table 3 together with 
the graviton bulk flight time and its corresponding brane time according to 
the equation (8). 

We have classified all cases according to the sign of the M parameter and 
to whether we are in dS or AdS bulks. Moreover, we studied the zero charge 
black hole as well as the Reissner-Nordstrom-type solutions, namely eight 
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cases. 



M 


k 


I 2 


Bulk 


DW 


MDB 


RDB 


Geodesic 


+ 


1 


+ 


AdS-Schwarzschild 


^r H /G 


-> r ff 


-> r ff 


- > th* / G 


+ 


0,-1 


+ 


AdS-top. black hole 


-> r H 


-> r# 


-> r ff 


-> rjj* 




0,1 


+ 


AdS-naked singularity 


G 


-> 


-> 






-1 


+ 


AdS-top. black hole 


-> r H 


-> r H 


-> r ff 


-> r ff 


+ 


1 




dS-Schwarzschild 


-> r H /r c 


-> rff/^c 


-> r.ff/r c 


-> r H /r c *^ 


+ 


0,-1 




dS-cosm. singularity 


no solution 


no solution 


no solution 


no solution 




0,±1 




dS-naked singularity 


— » r c 


^0/r c 


^0/r c 


-» 0/r c *t 



Table 1: Scale factor and geodesies evolution (uncharged case). The arrow 
indicates the behaviour tendency, which depends on the initial conditions. G 
means growing behaviour. The * or the f in the last column indicates the 
possibility of shortcuts for the matter and radiation-dominated branes or the 
domain wall, respectively. 




Figure 2: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M > 0, Q 2 = and I 2 > 0. 
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R(t) for M<0, Q 2 =0, k=l, 1 2 >0 R(t) for M<0, Q 2 =0, k=l, 1 2 >0 



(Zoom of the region of interest) 




Figure 3: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M < 0, Q 2 = and I 2 > 0. 



First, let us consider the zero charge case. The following results have 
also been summarized in table 1 for a better comprehension of the several 
solutions. 

4.1 M > 0, Q 2 = 0, I 2 > 

For positive both M and I 2 when k — 1, the bulk is anti de Sitter- 
Schwarzschild solution, while for fc ^ 1 we have a "topological" black hole 
with a flat or hyperbolic event horizon in an asymptotically anti de Sitter 
space. This case is shown in figure 1(a). 

For a domain wall (u = —1) and k = 1 the brane falls into the event hori- 
zon th when the initial condition for the brane position is near th, otherwise 
lZ(t) grows. On the other hand, for k ^ 1 the brane always falls into r#. 

In the case of a matter-dominated brane (u = 0) the brane falls into 
the event horizon for any k. This behaviour is also verified for a radiation- 
dominated brane (a> = 1/3). 

Let us see now the geodesic behaviour. When k — 1 and the initial 
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R(t) for M>0, Q 2 =0, k=l, 1 2 <0 fi(t) for M>0, Q 2 =0, k=l, 1 2 <0 




(112345 2 4 6 8 



t t 

Figure 4: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M > 0, Q 2 = and I 2 < 0. 



condition on the brane is chosen near rjj, the geodesies fall into the event 
horizon. Some of them follow the domain wall at the beginning and others 
leave the matter or radiation-dominated branes returning after a short time 
(shortcut geodesies) and leaving again the branes to fall into the singularity. 
Conversely, if the initial condition is taken far from the event horizon, the 
geodesies grow accompanying the domain wall initially, though their growth 
is slower; however, they never return to it after their decoupling. When 
k 1, independently of the initial condition the geodesies behave in the 
same way as the previous case when the initial condition is taken near r# 
and we again find some shortcuts for the matter and radiation-dominated 
branes before all the geodesies fall into the event horizon. Some results are 
shown in figure 2. 

4.2 M < 0, Q 2 = 0, I 2 > 

As we can see from figure 1(b), for negative M and positive I 2 when k = — 1, 
h describes a topological black hole in an asymptotically anti de Sitter space 
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«(t) for M<0, Q 2 =0, k=l, 1 2 <0 



fi(t) for M<0, Q 2 =0, k=l, 1 2 <0 




Figure 5: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M < 0, Q 2 = and I 2 < 0. 



with hyperbolic event horizon. If k — 1, there is a timelike naked singularity 
and the metric is asymptotically anti de Sitter. 

The solutions of the brane equation of motion (28) display the following 
features. In the case of a domain wall, when k = —1, the brane falls into the 
event horizon; while if k ^ —1, the solution for 7Z(t) grows. 

For (u = 0) the matter-dominated brane falls into the naked singularity 
when k ^ — 1 and into the event horizon if k — — 1. The radiation-dominated 
brane displays the same behaviour. 

On the other hand, the solutions of equation (29) show that when k ^ —1, 
the geodesies grow slower than the domain wall. Some of them are shortcuts 
for matter or radiation-dominated branes since they leave and return to them 
before the branes reach the naked singularity. Furthermore, if the geodesic 
(negative) initial velocity is big enough, it can fall into the naked singularity. 
Besides, when k = —1, all the geodesies fall into the event horizon. We can 
see some results in figure 3. 
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dS-naked singularity 
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-» r c 


^0/r c 


- r c *t/0 
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dS-Reissner- Nordstrom 


— ► r H /r c 


— > r_fj-/r c 


— > r ff /r c 


— > r H * /r c *^ 
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0,±1 




dS-naked singularity 


—> r c 


^0/r c 


^0/r c 
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Table 2: Scale factor and geodesies evolution (charged case). B and G mean 
bouncing and growing behaviour, respectively. The * or the f in the last 
column indicates the possibility of shortcuts for the matter and radiation- 
dominated branes or the domain wall, respectively. 



4.3 M > 0, Q 2 = 0, I 2 < 

In the case of positive M, negative I 2 and k — 1, the metric is de Sitter- 
Schwarzschild with event and cosmological horizons given by the zeros of 
h(7t). If k ^ 1, there is a cosmological singularity at 1Z = in the asymp- 
totically de Sitter background. See figure 1(c). 

In the case of a domain wall for k — 1, the solutions for lZ(t) converge 
either to the event or the cosmological horizon. For k ^ 1, there is no solution 
since 1Z turns to be a time coordinate. The same behaviour is observed for 
matter and radiation-dominated branes. 

Since the solutions to the brane equation of motion just appear when 
k — 1, the solutions of the geodesic equation only have physical meaning in 
this case. In this way we found that the geodesies reach either the event or 
the cosmological horizon. In the latter case there are many shortcuts for the 
matter or radiation-dominated branes as well as for the domain wall before 
all of them reach the cosmological horizon. Some results are illustrated in 
figure 4. 
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«(t) for M>0, QVo, k=0, 1 2 >0 



R(t) for M>0, QVo, k=0, 1 2 >0 




Figure 6: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M > 0, Q 2 7^ and I 2 > 0. 



4.4 M < 0, Q 2 = 0, I 2 < 

When M and / 2 are both negative, the bulk metric displays a timelike naked 
singularity in asymptotically de Sitter space. There is also a cosmological 
horizon which geometry is determined by k as we can see from figure 1(d). 

In the case of a domain wall for any k the solution TZ(t) converges to the 
cosmological horizon. 

When (u = 0), the matter-dominated brane either falls into the naked 
singularity or converges to the cosmological horizon independently of the 
value of k. This same result was found for a radiation-dominated brane. 

On the other hand, the geodesies either fall into the naked singularity 
when the initial velocity is negative enough or converge to the cosmological 
horizon for any k. In the latter case we found several shortcuts for the domain 
wall and the branes. We can see some results in figure 5. 

Now let us consider the charged solutions. As in the uncharged case, we 
also show a summary of our results in table 2. 
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Figure 7: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M < 0, Q 2 ^ and I 2 > 0. 



4.5 M > 0, Q 2 ^ 0, I 2 > 

When M and Z 2 are both positive, for k — —1 the metric describes a topolog- 
ical charged black hole in asymptotically anti de Sitter space with hyperbolic 
horizon. While for k ^ — 1, there is a timelike naked singularity and the bulk 
is asymptotically anti de Sitter (see figure 1(e)). Note that if the charge is 
small, all the metrics describe topological charged black holes in anti de Sitter 
bulks. In particular, when k — 1, we have anti de Sitter- Reissner- Nordstrom 
bulk. 

Let us see the solutions of the brane equation of motion. In the case 
of a domain wall there are different behaviours for each value of k. When 
k = 0, 7Z(t) bounces for the initial condition 71(0) greater than or equal 
to the minimum of h(7Z); otherwise, it diverges to infinity. When k = 1, 
7Z(t) diverges to infinity. When k — — 1, the domain wall falls into the event 
horizon. For small charges 7Z{t) converges to the event horizon except when 
k — 1, where it can also grow or diverge to infinity after a certain initial 
value. In the limit Q — > we also verified that our solutions converge to 
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R(t) for M>0, QVo, k=l, 1 2 <0 R(t) for M>0, Q 2 «l, k=0, 1 2 <0 



(Zoom of the region of interest) 




Figure 8: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M > 0, Q 2 7^ and I 2 < 0. 



those in section 4.1. 

For a matter-dominated brane {u = 0) we again found different be- 
haviours for each k. When k = 0, 7Z(t) bounces for the initial condition 
7Z(0) less than or equal to the minimum of h(TZ); otherwise, it falls into the 
naked singularity. When k — 1, 7Z{t) goes to the naked singularity. When 
k = —1, it falls into the event horizon. For small charges the brane always 
falls into the event horizon. Again when Q — > we recover the solutions in 
section 4.1. 

In the case of a radiation-dominated brane (u = 1/3) the behaviour of 
the solution is somewhat similar to the matter-dominated brane case. When 
k = 0, 7Z(t) bounces for the initial condition 71(0) less than or equal to the 
minimum of h(7Z), but it falls into the naked singularity when the initial 
condition is near or when it is greater than the minimum of h(lZ). When 
k = 1, 7Z(t) goes to the naked singularity. When k = —1, it falls into the 
event horizon as the previous case. For small charges we have the same 
behaviour as in the matter-dominated brane case. 
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R(t) for M<0, Q 2 #0, k=l, 1 2 <0 R(t) for M<0, Q^O, k=l, 1 2 <0 
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t t 



Figure 9: Scale factor evolution for domain wall, matter and radiation- 
dominated branes and geodesies when M < 0, Q 2 ^ and I 2 < 0. 



Now let us investigate the geodesic behaviour. When k = 0, the geodesies 
bounce for the initial condition less than or equal to the minimum of h(JZ) 
in the same way as matter or radiation-dominated branes producing several 
shortcuts in all brane cases. However, when the initial condition is greater 
than the minimum of h(lZ), the geodesies can continue to bounce or escape 
after one oscillation or even fall into the naked singularity as a result of 
increasing the initial (negative) velocity; in this case we just found shortcuts 
for the domain wall as long as the geodesies complete at least one oscillation. 
For small charges an event horizon is formed and the geodesies fall into it but 
there are still some small shortcuts for matter or radiation-dominated branes. 
When k — 1, the geodesies grow slower than the domain wall producing 
shortcuts just for matter or radiation-dominated branes. For small charges 
the geodesies either diverge or fall into the event horizon. When k — — 1, 
all the geodesies fall into the event horizon as the branes and there are no 
shortcuts. For small charges the geodesies fall into the event horizon but 
can yield some shortcuts in the very beginning of their paths for radiation 
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Table 3: Bulk time t, brane time r, time delays At and ratio between graviton 
and photon horizons (7/7 for shortcut geodesies. 



or matter-dominated branes. Some of these cases appear in figure 6. 
4.6 M < 0, Q 2 ^ 0, I 2 > 

For negative M, positive I 2 and k — — 1, the metric describes a topological 
charged black hole in asymptotically anti de Sitter spacetime with hyperbolic 
horizon. While for k 7^ — 1 there is a timelike naked singularity and the bulk 
is asymptotically anti de Sitter as is shown in figure 1(f). 
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In the case of u — — 1 and for k — — 1 the domain wall falls into the event 
horizon. If in turn k ^ — 1, the solution 1Z(t) diverges to infinity. 

When we consider a matter-dominated brane for k = —1, the brane 
behaves in the same way as a domain wall; however, if k ^ — 1, it falls 
into the naked singularity. This same behaviour is observed for a radiation- 
dominated brane. 

As for the solutions to the geodesic equation, when k = — 1, the geodesies 
fall into the event horizon and some of them can produce shortcuts for radia- 
tion or matter-dominated branes before their falling. If k ^ — 1, the geodesies 
grow and yield shortcuts for matter or radiation-dominated branes. Some re- 
sults are shown in figure 7. 

4.7 M > 0, Q 2 ^ 0, I 2 < 

In the case of positive M and negative I 2 the metric describes a timelike naked 
singularity in asymptotically de Sitter space with a cosmological horizon as 
we see in figure 1(g). However, for a small charge the k = 1 metric turns 
out to be de Sitter-Reissner-Nordstrom with Cauchy, event and cosmological 
horizons as it is shown in figure 1(h). 

For a domain wall the solution TZ(t) converges to the cosmological horizon 
for any k. For a small charge this behaviour also applies except for k — 1, 
when the domain wall can also fall into the event horizon. 

The case of a matter-dominated brane displays an interesting behaviour. 
When k — 1, 7Z(t) bounces when the initial condition 71(0) is less than or 
equal to the only saddle point in h(TZ); otherwise, it falls into the naked 
singularity or converges to the cosmological horizon if the initial condition 
7^.(0) is very near it. When k ^ 1, the solutions converge to the cosmological 
horizon. For a small charge the solutions converge to the cosmological horizon 
and also to the event horizon in the case k = 1. 

When we consider a radiation-dominated brane, we see that it either falls 
into the naked singularity or converges to the cosmological horizon for any 
k. For a small charge the solutions behave in the same way except if k — 1 
when the brane converges either to the cosmological or the event horizon. 

Let us see the geodesic behaviour. For any k all the geodesies converge to 
the cosmological horizon and produce several shortcuts for all brane cases. 
When the charge is very small, an event horizon appears in the case k = 1 
and the geodesies can either converge to the cosmological or the event hori- 
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zons. The geodesies converging to the cosmological horizon produce shortcuts 
for all brane cases, while the geodesies falling into the event horizon yield 
shortcuts just for matter or radiation-dominated branes. On the other hand, 
when k ^ 1, the geodesies converge to the cosmological horizon after produc- 
ing some shortcuts for all brane cases, unless their initial (negative) velocity 
reaches a threshold after which they fall into the naked singularity. Some 
results are illustrated in figure 8. 

4.8 M < 0, QV 0, I 2 < 

When M and I 2 are both negative, the metric describes a timelike naked 
singularity in asymptotically de Sitter space for any k (see figure l(i)). 

As for the solutions of the brane equation of motion, in the case of a 
domain wall all the solutions converge to the cosmological horizon. 

When we consider a matter-dominated brane, TZ(t) either falls into the 
naked singularity or converges to the cosmological horizon for any k. The 
same behaviour applies for a radiation-dominated brane. 

As for the geodesies, they converge to the cosmological horizon and pro- 
duce several shortcuts for all brane cases for any k. However, when the 
initial condition is taken near the singularity and the initial velocity is nega- 
tive enough, the geodesies can reach the singularity and no shortcut appears. 
We can see some results in figure 9. 

5 Discussion and Conclusions 

In the present work we have studied the behaviour of a brane embedded in a 
six-dimensional de Sitter or anti de Sitter spacetime containing a singularity 
covered by at least one horizon in the case of black hole type solutions or a 
timelike naked singularity. The system of equations describing this behaviour 
from the point of view of an observer in the bulk appears to be highly nonlin- 
ear. Before numerically solving this system we have considered fluctuations 
about a fixed brane position in order to have a better insight of the whole 
problem. We have concluded that the case of a flat domain wall with van- 
ishing effective cosmological constant reproduces the equation of motion for 
a free scalar field as in five dimensions [15, 18]. The "equilibrium" position 
can be chosen arbitrarily but there is no stability. 



23 



By solving the full nonlinear system we found different behaviours for the 
several scenarios appearing due to all the combinations of M, Q 2 , k and I 2 
taken into account. We chose some typical values for cu, i.e. domain wall, 
matter and radiation-dominated branes, in order to illustrate the solutions. 
The results show branes getting away from the singularity, falling into it, con- 
verging to cosmological horizons when they exist or even bouncing between 
a minimum and maximum values. 

The bouncing behaviour found in some of our solutions of the brane equa- 
tion of motion appears to be in good agreement with the recent investigations 
in five dimensions [19], where universes bouncing from a contracting to an 
expanding phase without encountering past and/or future singularities ap- 
pear. In this way these results could provide support for a singularity-free 
cosmology or to the so-called cyclic universe scenarios [20]. 

Finally, we also studied the geodesic behaviour in every scenario found in 
the present work. Contrarily to the case of a static brane, where shortcuts 
appeared under very restrictive conditions [6] , the present model of a dynamic 
brane embedded in a static bulk displays shortcuts in almost all cases and 
under very mild conditions. Moreover, despite the fact that the time delay 
between graviton and photon flight time is not percentually so big as in 
other models [17] (what is also evident from the ratio between graviton and 
photon horizons), it exists and can eventually be measured by the brane 
observer, although further considerations are certainly needed in a stricter 
realistic model. On the other hand, the fact that shortcuts are abundant in 
the studied setups lends further support to the idea of solving the horizon 
problem via thermalization by graviton exchange [17, 13]; however, we should 
stress that this is not a proof of the solution of the problem yet. 
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